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Abstract 

Two-surface wave decay (TSWD), i.e. the parametric excitation of electron surface waves, was 
recently proposed as an absorption mechanism in the interaction of ultrashort, intense laser pulses 
with solid targets. We present an extension of the fluid theory of TSWD to a warm plasma which 
treats boundary effects consistently. We also present test-particle simulations showing localized 
enhancement of electron acceleration by TSWD fields; this effect leads to a modulation of the 
current density entering into the target and may seed current filamentation instabilities. 

PACS numbers: 52.38.-r: 52.38.Dx 
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I. INTRODUCTION 



The excitation of electron surface waves (ESWs) is a possible route to collisionless energy 
absorption in plasmas produced by the interaction of sub-picosecond, high intensity (typi- 
cally > 10 16 W/cm 2 ) laser pulses with solid targets. This regime is relevant for applications 
such as generation of ultrashort X-ray pulses, either as uncoherent thermal emission and 
references therein) or high laser harmonics 2], or production of energetic electrons and ions 

□ □□□ 

Linear mode conversion of the laser wave into an ESW is not possible at a planar plasma- 
vacuum interface, because the phase matching between the ESW and the incident wave is 
not allowed; the process can take place for specially tailored density profiles, e.g. grating 
targets with a density profile modulated at the surface with a wavevector k p . In this case, 
the condition for the excitation of a ESW with wavevector k s (parallel to the plasma surface) 
and the same frequency ujl of the laser pulse is 

k L sm8 = k s + k p , (1) 

where k^ = Ul/c, 9 is the angle of incidence of the laser pulse. Experimental investigations 



a. 



of laser absorption and X-ray emission in grating targets are described, e.g., in Refs.[ 

Nonlinear mode conversion, i.e. parametric excitation of two ESWs is possible also in 
a simple "step" profile. We name this process "two-surface wave decay" (TSWD) 
Examples of similar TSWD processes were previously considered in regimes other than 
intense laser-plasma interactions and in the electrostatic limit only 0, Q|. The general 
phase matching conditions for such a three-wave process are given by 

k = k + + k_ , luq = uj + + co>_, (2) 

where k , u and k±, u± are the wavevector and frequency of the "pump" wave and the two 
ESWs, respectively. If the electric field of the laser pulse acts as a "pump" for TSWD, then 
uo = u>l, ko = kisinO, and equation (J2J) implies that two sub-harmonic ESW are generated 
with frequency around u/2 ( u u — > uj/2 + cj/2" TSWD). However, also the v x B term of 
;he Lorentz force may excite TSWD. This was observed in particle- in-cell (PIC) simulations 



111 ] for normal laser incidence; in this case, the v x B force drives ID oscillations at the 
plasma surface with frequency 2u>l] after a few laser cycles the overlap of a standing wave 
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with frequency ujl was observed. This is a clear signature of a u 2u — > oo + uj n process leading 
to two counterpropagating ESWs both having the frequency of the laser pulse. 

According to the theoretical model jsl Q], the maximum growth rate of the 2uo — > uo + uo 
case is found for normal incidence. As said above, TSWD does not need a structured 
target; however, it is worth noticing that in a grating target the wavevector k p of the surface 
modulation [eq.flT}] equals the wavevector of the two ESWs excited by the 2uo — * uo+uo process 
at 9 = 0. Hence, TSWD is enhanced in such a grating target. It is then interesting to notice 
that PIC simulations reported in Ref. Q], showing the generation of two countepropagating 
ESWs in grating targets at 10 16 W/cm 2 , can now be interpreted as an evidence of TSWD 
seeded by the surface grating. 

At very high intensities (> 10 18 W/cm 2 ), simulations |h| show that TSWD enters a 
strongly nonlinear regime leading to strong rippling of the plasma surface. This might be 
relevant to the surface instabilities which have been observed in experiments [131 and play a 

_ __ — Eflua 

appear to grow even for pulses of few tens of fs [Iq and thus must be of electronic nature, 
since ions do not move on such a fast scale. Simulations also show that nonlinear TSWD 
affects fast electron generation. This effect is investigated in section IIHI by means of test- 



particle simulations. Before that, in section |H] we present an improvement of the model of 
Ref. s| where temperature effects have been included. In both cases, we restrict ourselves to 
the case of the 2uj — > uo + uo TSWD at normal laser incidence, and we neglect collisions for 
simplicity and because collective processes dominate absorption at intensities > 10 16 W/cm 2 . 



II. TEMPERATURE AND BOUNDARY EFFECTS ON TSWD 

In Ref. p| the TSWD growth rate was calculated for a step-like density profile [rij = 
rioO(x), being Q(x) the Heaviside step function] using an Eulerian, fluid model with immobile 
ions and using the cold plasma approximation. We adopted the following expansion for all 
fields 

f(x,y,t) = f t ( x ,t-ysm6/c) + ef (x)e^ t -y si ^ + 

e 2 [f + (x)e ih+y - iuJ+t + f.{x)e ik - y - ju -% (3) 
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where e is a small expansion parameter, and / stands for either the electron density or 
velocity or for the EM fields in the (x, y) plane. The first term (fi) of eq.Q includes zero- 
order, unperturbed fields or oscillating fields that are non-resonant with the excited modes; 
the term f represents the "pump" field at the frequency ujq; the last term is the sum of 
two counterpropagating surface modes. For the 2u — > u> + u> process at 9 = 0, u = 2uj l 
and k + = —k-, and the zero- and first-order fields do not depend on y (i.e. they are "ID" 
fields). The coupling between the pump and the surface modes (of order e 3 ) originates from 
the nonlinear terms in the Euler equation (— en e v x B, — m e n e vVv) and the current density 
J = — en e v. 

The calculation for a warm, isothermal plasma proceeds very similarly 17[ and thus 
only the differences from the "cold" case and their consequences are discussed below, while 
details of the calculation will be published elsewhere. The only difference in the starting 
Maxwell-Euler systems of equations comes from the pressure term in the Euler equation for 
electrons: 

m e dv/dt = -en e (E + v x B) - k B T e Vn e , (4) 

where n e , v and T e are the density, fluid velocity and temperature of electrons. 

Assuming reflective conditions at the plasma boundary (x = 0), the charge density p = 
e(rii — n e ) ^ only in a surface layer with a thickness on the order of the Debye length. 
At T e = this corresponds to a surface charge layer [i.e. p ~ S(x), where 8(x) is the Dirac 
delta function], while the longitudinal velocity (v x ) is discontinous at the surface. This is 
relevant because in the T e = case nonlinear coupling terms involving the product of fields 
that are singular at x = occur. The pressure term removes such singularities and enables 
to evaluate all such terms correctly. Therefore, by calculating the TSWD growth rate and 
then taking the T e — > limit one achieves an improvement of the result obtained in the cold 
plasma case. 

The final result of the calculation is the 2a; — »• u> + u growth rate T shown in FigQ (the 



complete analytical expression is very lengthy [12 1 and is not reported here). With respect 
to the T e = case, the two divergences of 7 are both quenched for growing T e , but for 
different reasons. The divergence for oj — > uj p / y/2 is quenched since in this limit the ESWs 
have shorter wavelengths and are thus more affected by the thermal pressure that inhibits 
the formation of small-scale structures. The "pump" resonance at uj ~ uj p /2, due to the 
excitation of plasmons by the v x B force, is quenched by energy transport out of the skin 



layer because for T e ^ the resonant plasmon propagates into the overdense plasma. 



III. ELECTRON ACCELERATION 

At high intensity of the laser pulse, it is well known that most of the absorbed energy goes 
into "fast" electrons injected into the overdense plasma region. At normal incidence, due 
to the leading action of the v x B force, fast electron bunches are produced twice per laser 
cycle. The PIC simulations of Ref. showed that after the growth of nonlinear TSWD 
the generation of fast electrons was enhanced near spatial maxima of the standing surface 
oscillation (Fig|21 top). Thus, this effect may give an "imprint" for the formation of electron 
filaments, whose size and spacing would be close to (and scaling with) the laser wavelength 

nn 

as observed in other simulations |18l . [19J , and affect energy transport by fast electrons into 
the plasma. 

To investigate this effect further, we performed "test particle" simulations of electron 
motion in the overlapping "pump" and TSWD fields. In other words, we solve the equation 
of motion for electrons into a force field of the form given by eq.Q, i.e. the sum of an 
one-dimensional (ID) force of frequency 2uj (whose analytical expression is obtained from 
the theoretical model) and the force field from a standing 2D surface wave of frequency 
uj. A similar study, focused on the acceleration of electrons by a single ESW, was recently 
reported j^J. 

The ID and 2D fields vary in time as cos2co>£ and sin tut, so that the temporal maxima 
and minima of the ESW field are always coincident in time with maxima of the ID field, 



as it was found in theory and simulation ja, In what follows, x is the direction normal 
to the plasma surface located at x = (the plasma occupies the x > region) and y is the 
direction of propagation of ESW. The field amplitudes are chosen so that the system remains 
far from relativistic conditions and the ESW field can be considered as a perturbation. The 
ratio between the amplitudes of standing SW and the 2a; pump was varied as a parameter. 
For the simulations reported in this paper, in terms of normalized amplitudes, a = eE/m e cu> 
where E is the amplitude of the electric field, the values are = 0.2 for the laser field and 
a esw — 0.019 for the ESW field. The wavevectors of the laser wave and of the ESW are 
given by the known expressions for a cold plasma with density n e /n c = oj 2 /uj 2 = 5. To 
initialize the simulations, we gave the particles an initial position a;(0) > such that the 
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evanescent fields are negligible at that point, and an initial velocity 1^(0) < (on average 
v x = —0.1c) so that the particles move towards the surface region, where the ESW fields are 
localized. The particles are distributed in y uniformly over a region of width A s , where X s 
is the ESW wavelength. 

The (y,p x ) projection of the test particle phase space (FigJ21 bottom) looks similar to 
the one from the fully self-consistent PIC simulation (Fig0 top), showing that electron 
acceleration is indeed enhanced near spatial maxima of the standing wave. The overlap of 
the SW fields with the v x B-driven fields (homogeneous along the ^-direction) thus leads 
to a modulation along y of the longitudinal momenta. 

The enhancement (or quenching) of the accelerating field occurs once per laser cycle 
(because the ESW frequency equals the laser frequency), but with opposite phase between 
contiguous maxima of the standing ESW. This is shown by Fig|3] where the complete (x,p x ) 
phase space projection, showing electron bunches generated at 2ui rate, is compared with the 
same plot but including now only the electrons whose starting position lies within an interval 
of A/4 width around the y = A/4 spatial maximum of the ESW. It is thus evident how the 
most energetic electrons in the bunches penetrating in the x > region are generated near 
the ESW maximum and at u rate, i.e. with the ESW frequency. Taking only the electrons 
around the y = 3A/4 spatial maximum gives an almost identical picture, except that the 
most energetic bunches are now out of phase by an angle of 2tc/uj with respect to those 
coming from y = A/4. 

Although the amplitude of ESW fields is 0.1 times the pump field amplitude, the relative 
modulation of the longitudinal momentum p x is about 30%. This can be qualitatively 
explained at follows. The electrons acquire energy from an evanescent, oscillating field if 
their transit time across the region where the field has non-vanishing amplitude is shorter 
than an oscillation period; this is the condition for non-adiabaticity of electron motion 
with respect to the field. The ratio between the transit time and the oscillation period is 
approximately given by the parameter rj = L/vqT where L is the evanescence length, vq is 
the average velocity of the electron and T is the oscillation period. For the ESW, T = 2n/u) 
and Iesw — (c/uS)\/ot. — 2/ (a -^1) (where a = ujp/uj 2 ) while for the ID field at T — ir/oj and 



hui = (c/2u;)(l/v / tt — 1) and |8[. Thus, rjEsw/vzu = \J ( a ~ 2)/(a — 1) < 1 which means 
that the electron motion (for a given vq) is more non-adiabatic with respect to the ESW 
field rather than to the 2u oscillation. Thus, the contribution of the ESW in accelerating 
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(or decelerating) electrons is enhanced by the relatively low frequency and short evanescence 
length. 

The "imprint" effect of the standing ESW is also noticeable in the contourplot of the fluid 
velocity of electrons, shown in Fig0] The velocity field has been computed by averaging the 
velocity over a spatial grid, as in PIC codes, and over time. As a consequence the electric 
current in the overdense plasma region has a spatial transverse modulation with the same 
wavelength of ESWs. The simulations of Ref . [t| for "grating" targets also show a modulation 
of the fast electron current correlated with same period of the grating, that may be due to 
the local enhancement of the lon gitu dinal field by TSWD or to the geometrical "funnel" 



effect of the surface deformation 



2l|. Our test particle results indeed show that due to 



TSWD a modulated fast electron current can be produced even by a plane-wave pulse on a 
flat surface. 

From figures El and H3 ft is also found that only near ESW spatial maxima a few electrons 
are ejected into vacuum, i.e. in the x < region (in FigOl for these electrons the oscillation 
of p x vs. x is due to the effect of the v x B force in vacuum). Their origin is likely to 
be due to the l ong itudinal field component directed into vacuum that is associated to the 
ESWs. In Ref. |20| . the features of electrons accelerated in vacuum by the field of a single 
ESW are investigated. Note, however, that electrostatic back-holding fields are not self- 
consistently included in test particle simulations and thus the number of electrons escaping 
in vacuum is likely to be overestimated. Nevertheless, it is interesting to notice that in the 



PIC simulations at very high intensities III, ll2| "plumes" of electrons extending into the 



vacuum region are generated near the maxima of the standing ESW. A "plume" structure 
in the vacuum region is also evident in the velocity field shown in Fig0] 
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FIG. 1: Growth rate of the 2lo — > oj + uj process vs. normalized electron density n e /n c = Up/ui^ 
(x-axis) and temperature \/kBT e /m e c 2 (labels). The dashed line is the "cold" result previously 
obtained [8(]. The growth rate is normalized to a 2 L uj where = eEi/m e ajJL is the dimensionless 
field amplitude of the laser pulse. 

FIG. 2: Top: (y,p x ) phase space projections from PIC simulations at two subsequent times, 
showing electron acceleration localized near spatial maxima of the growing surface oscillations. 
Bottom: (y,p x ) phase space projection from test particle simulations, showing similar features. 
The "black stripe" around v x = —0.1 represents electrons that have not been reached yet the 
surface region at the time shown. 

FIG. 4: Contours of time-averaged fluid velocity from test particle simulations, showing a transverse 
modulation due to the enhanced electron acceleration near ESW maxima (y ~ A s /4, y ~ 3A s /4). 



FIG. 3: Left: phase space (p x ,x) projection (integrated over y) from test particle simulations, 
showing fast electron jets penetrating into the plasma region (x > 0) where they move ballistically 
(p x oc x). The jets are produced twice per laser cycle. All particles in the simulation are included 
in the plot. Right: same as left side, but restricted to particles around y = A s /4 (spatial maximum 
of the standing surface wave), showing enhanced acceleration around this position once per laser 
cycle. In both figures electrons propagating into vacuum (x < 0) are also evident. They are found 
only near ESW maxima {y ~ A s /4, y ~ 3A s /4). 
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